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The main idea of this article is based on my previous articles (references [1], [2], [3]). In this work by 
introducing a new mathematical approach based on the algebraic structure of integers (the domain of 



a 

D integers), and assuming the "discreteness" of physical quantities such as the components of the relativistic 

W) 

• n-momentum, we derive all the mathematical laws governing the fundamental forces of nature. These 

obtained laws that are unique, distinct and in the form of the complex tensor equations, represent the force 

• i-H 

of gravity, the electromagnetic force, and the (strong) nuclear force (and only these three kinds of forces, 
for all dimensions D >2). Each tensor equation contains the term of the mass mg (as the invariant mass of 

rv 

i the supposed force carrier particle), as well as the term of the external current (as the external source of 

the force field). In some special cases, these tensor equations are turned into the wave equations that are 
similar to the Pauli and Dirac equations. In fact, the mathematical laws obtained in this paper, are the 
corrected and generalized forms of the current field equations including Maxwell equations, Yang-Mils 

•l-H 

equations and Einstein equations, as well as (in some special conditions) Pauli equation, Dirac equation, 



and so on. A direct proof of the absence of magnetic monopoles in nature is one of the outcomes of this 
research, according to the unique formulations of the laws of the fundamental forces that we have derived. 



1. Introduction 

One of the biggest ontological questions might be that: "Why the universe and the fundamental forces that 
are acting in it, are in the way and form, which we realize them?"; the forces that are the causers and the 
governors of all of the movements, the actions and reactions in the physical world. In this article, we are 
going to consider this question by a logical approach. 

This paper is based on my previous articles ([1], [2], [3]), as well as my thesis work, (1997) [4], (but in a 
new expanded framework). Here by representing a new axiomatic mathematical approach, as well as 
assuming the discreteness of the components of the relativistic n-momentum, the basic laws governing 
the fundamental forces of nature will be derived. These laws are unique and in the form of tensor 
equations. The main results of this paper include: 
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1-1. The mathematical laws governing the fundamental forces of nature (for all dimensions): 



*VW=o, (i-i) 



D [A Z MV]p =0, (2-1) 



*VV]=0, (3-D 
D;F»=-J { v e) (3-2) 



where F is the field tensor of the electromagnetic force, Z is the field tensor of the nuclear (strong) 
force, and R M vpa is the Riemann tensor of the gravitational field, and 

*>,=v„+^*„. ^;=v,-^^ (4) 



// = 0: 



V ^ (5) 
M*0: k M =0, 

where mo is the invariant mass of the (supposed) force carrier particle. 

As we will show in the section 3., equations (1-1) - (1-2) (representing the gravitational fields), also could 
be written as follows, 

for two dimensional case (D = 2) we have 



R-A = -SxT (6) 
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R MV = -y 2 (Sfdr)8 MV + y2Ag MV (7) 

jffl 

and for higher dimensions (D > 2) we have 

ifti 

R MV -yiR8 MV = -^T MV -^K MV -A 8fJV (9) 

where K MV = VJc v . 

1-2. According to the unique structure of the tensor equations (3-1) and (3-2) that in fact, are the 
corrected and generalized form of the Maxwell equations, we will conclude that there are not any 
magnetic monopole in nature. 



In the next section, we describe the bases of the linearization theory, and in the section 3. we show its 
applications in physics. 



2. Theory of Linearization in the Domain of Integers: As a New Axiomatic 
Mathematical Approach 

The algebraic axioms of the domain of integers Z , with binary operations (+,x), usually are defined as 
follows [5]: 

-Closer: a k +a l sZ, a k xa,sZ (10) 

-Associativity: a k + (a, + a r ) = (a k + a t ) + a r , a k x(a l xa r ) = (a k xa,)xa r (11) 

- Commutativity: a k +a l =a l +a k , a k xa t = a t xa k (12) 

- Existence of an identity element: a k + 0 = a x , a k x 1 = a k (13) 
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- Existence of inverse element (for addition): a k + (-a k ) = 0 (14) 

- Distributivity : a k x (a l +a r ) = (a k xa l ) + (a k xa r ), 

(a k +a,)xa r =(a k xa r ) + (a,xa r ) (15) 

- No zero divisors: (a k = 0va, = 0)oa k xa, = 0, (16-1) 

equivalently, the axiom (16-1) could be defined as: 

[(a,xm, =0,mj ^ 0) v (a 2 xm 2 = 0,/?% ^0)v...v 

v (a r x m r = 0, m r * 0)] oa 1 xa 2 xa 3 x...xa r = 0. (16-2) 

If we just suppose [a] lxl (= a),[b] lxl (= b),[c] lxl (= c),... e Z M (= Z), then equivalently, the axioms 
(10) - (15) could also be written by quadratic matrices (with integer components) as follows: 

-M k =[m ki ], m,.eZ, 3n e N :/, ; = 1,2,3,..., n , M 1 ,M 2 ,M 3 ,... e Z„ xn , 

- Closer: M, +M, e Z nx „, M k xM t G Z nxn (17) 
-Associativity: M 4 + (M, +M r ) = (Af t +M Z ) + M r , M k x(M,xM r ) = (M k xM l )xM r (18) 

- Commutativity (for addition) : M k +M l =M l +M k (19-1) 

- Property of the transpose for matrix multiplication: 

(M k xMf = M, T xM k ( 19-2) 

T 

where M k is the transpose of matrix M k . 

- Existence of an identity element: M k + 0 = M k , M k x I nxn = M k (20) 

- Existence of the inverse element (for addition): M k + (—M k ) = 0 (21) 

- Distributivity: M k x(M,+M r ) = (M k xM,) + (M k xMJ, 

(M k +M l )xM r =(M k xM r ) + {M l xM r ); (22) 

From the axioms (10) - (15), we can obtain the axioms (17) - (22) and vice versa. 

In this article, we introduce the following algebraic axiom as a new property of integers, and we add it to 
the axioms (17) - (22) (this new axiom is someway the generalized form of the axiom (16-2), and in fact, 
the axiom (16-2) will be replaced with the axiom (23)): 
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"' If we assume the algebraic form F(b pq ) = / i \ \ b pq and the nxn quadratic matrices A k = [a k ] , 

q=l p=\ 

where H, „„ are some coefficients, and 

s r 

a k„ = HH H k,, Pq b p^ b pq ,H kiiPq ^Z (=Z ul ), i,j = 1,2,3,..., n , k =1,2,3,..., r , 
p = 1,2,3,..., r , q = 1,2,3,. ..,s , 
then we have the following axiom: 

3M k e Z nxn , [[(AjXMj = 0,Mj * 0) v {A 2 xM 2 = 0,M 2 * 0) v ... v xM r = 0,M r * 0)] a 
a (A, x A 2 x A 3 x ...x A, = F(^)/„ X J] F(b pq ) = 0 ."' (23) 

Remark 1. In (23), according to the arbitrariness of all of the components of the nxn matrix M k , 

without loss of generality, we may replace thenxn matrix M k with a nxl matrix M k , in each of the 

equations A k M k = 0 (with the same condition M k ^ 0, but only with the number of arbitrary 
components) . 

Note that the integer elements are the "linear" forms of the integer elements b pq . 

We can obtain the axiom (16-1) (or its equivalent, the axiom (16-2)) from the axiom (23), but not vice 
versa. Only for the special case n — l, the set of axioms (17) - (23) becomes equivalent to the set of 
axioms (10) - (16-2). Definitely, the axiom (23) is a new axiom and in this section and the next section 
we'll demonstrate some of its outcomes and applications. 

Generally, there are standard and specific methods, approaches and procedures for considering and 
solving the linear equations in the set of integers [7]. Since the necessary and sufficient condition for an 
equation of the r th order such as F{b pq ) = 0 (in the domain of integers) is the transforming (or 

converting or in fact, "linearizing") it into a system of linear equations of the type A k M k =0 

(according to the axiom (23), where M k ^ 0, M k : nxl matrix), naturally, the main application of the 

axiom (23) will be the transforming the higher order equations into the corresponding systems of linear 
equations. In this section, based on (23), essentially, we'll obtain the systems of linear equations that 
correspond to the second order equation of the form: F(b pq ) = 0, and also some of the higher order 
equations. In the methodological point of view, firstly, for obtaining and specifying a system of linear 
equations that corresponds to a given equation of the type F(b pq ) = 0 (defined in (23)), we assume and 

consider the minimum value for tl (the size number of nxn matrices A k ). Secondly, by replacing the 
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components of the matrices A k with the linear forms a k = ^^ H lnn b na , we calculate the product 



q=l p=l 



J~\_\ , and then we put it equal to the matrix F(b pq )I nxn . Then using this (obtained) equation, we can 



calculate the coefficients H k pq (which are independent of elements b pq ). Through this, easily, the 
coefficients H k pq are calculated and obtained by routine and standard methods of solving the equations 
(in the set of integers). Thirdly, the algebraic forms 

=zn*„ cm) 

q=\ p=l 

via some certain rules and linear transformations, could be transformed into the algebraic forms of the 
type 

S r 

G(c v c 2 ,c 3 ,...,c s )= FK ( 25 ) 

!',,/ 2 ,i' 3 ,.../,=l p=l 

s r 

In continuation, by some examples we will show how the forms F(b na ) = YY\^pq cou ^ be 

q=l p=l 

transformed into the forms G(c 1 , c 2 , c 3 ,..., c s ) (trough some linear transformations). Furthermore, as 

we'll show that, exceptionally, the second order forms of (24) could be transformed into the following 
quadratic forms (by similar linear transformation) 

s 2 s 2 

G(c 1 ,c 2 ,c 3 ,...,c J ,d 1 ,d 2 ,d 3 ,...,d,) = £ ^11^ " 2 B hhIl d i P W 

i 1 ,/ 2 =l p=\ ^ ,?2=1 P=l 



Remark 2. Concerning the formula (24), it is easy to show that 

s r s r s r 

(Xnvw),=°' znwr + i) 9 = °) ^ nn b p q << c (r + i ) q+d (r+1)q )=o , 

q=l p=l q=l p=l q=l p=l 

(24-1) 

^ Ell V(-D, = 0 ~ ill b P Mr +l) q) = 0- (24-2) 

q=l p=l q=l p=l 



where the parameter t is an arbitrary integer, with condition t ^ 0 . Below we write the systems of linear 
equations that correspond with some special cases of the following equation (that according to the axiom 
(23), one system for each case is enough): 
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q=l p=\ 



(24-3) 



that has been indicated in (23). Some special cases of (24-3), which we will consider below, in particular 
include the second order equations with the different number of the elements, as well as some of the 
other higher order equations. 

For s = 1,2,3,... , r = 2, equation (24-3) becomes as follows, respectively, 

Zri^ = *>iAi = 0, (27) 



q=\ p=l 
2 2 



zn^ =& iAi+M>22=o> 

q=\ p=\ 

3 2 

Zn^P? =b uK + b lAl + &L3&23 = 0, 
q=\ p=l 

4 2 

Zn^P? = ^ 1*21 + *12^22 + ^23 + ^At = 0, 



5 2 



Zn*/-? = b l 1*21 + *12^22 + &L 3&23 + KK + ^25 = 0; 
q=\ p=l 



(28) 



(29) 



(30) 



(31) 



Now, a matrix equation (here we mean a system of linear equations) corresponding to (27) (according to 
(23)) is 



~e 0 0" 




_° /o. 


m 2 



-0 



where e Q = b n , f Q = b 2l ; 

Similarly, for (28) we have the following matrix equation 



0 


0 


?o 


/i 


m l 


0 


0 




/o 


m 2 


fo 


fl 


0 


0 


m 3 






0 


0 


m 4 



= 0 



(32) 



(33) 
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where e 0 =b n ,f 0 =b 2V e l =b u ,f l = b 22 ; 

Using (33) we may get equivalently, the following matrix equation for (28) 



e 0 /l 


m l 


= 0 




m 2 _ 





(34) 



where e 0 =b n ,f 0 =b 2V e 1 =b u ,f 1 = b 22 ; 

A system of linear equations corresponding to (29) is 



0 


0 


0 


0 


^0 


0 


-e 2 


fx 


m l 


0 


0 


0 


0 


0 




-?x 


-/ 2 


m 2 


0 


0 


0 


0 


h 


/l 


fo 


0 


m 3 


0 


0 


0 


0 






0 


/o 


m 4 


-f 0 


0 


-h 


-«1 


0 


0 


0 


0 


m 5 


0 


-fo 


fx 




0 


0 


0 


0 


m 6 


e 2 






0 


0 


0 


0 


0 


m 1 


/i 




0 




0 


0 


0 


0 


_m 8 



where e 0 — £>j p /" 0 — b 2V e l — b l2 , f — b 22 ,e 2 — b lv f 2 — Z? 23 ; 

from (35) we can obtain the following matrix equation for equation (29) 



e Q 0 -e 2 f 
0 e o ~ e x ~fl 

h A fo o 

_-e l e 2 0 f 0 
where e 0 = b x v f 0 = b 2V e l = b u , f = b 22 ,e 2 = b 13 ,f 2 = b 23 ; 

Similarly, how we obtained the matrix equations (35) and (36) for equations (28) and (29), the matrix 
equations corresponding to (30) and (31) are obtained as follows, 

for (30) we get 



m. 



m 0 



m. 



= 0 



(36) 
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e 0 


0 


0 


0 


0 


-e 3 


e 2 


fx' 


m l 


0 




0 


0 


^3 


0 


-«i 


h 


m 2 


0 


0 




0 




«i 


0 


h 


m 3 


0 


0 


0 


^0 


-fx 


-/ 2 


~f 3 


0 


m 4 


0 


-U 


h 


e l 


fo 


0 


0 


0 


m 5 


/a 


0 


-fx 




0 


fo 


0 


0 


m 6 


-h 


/i 


0 


e 3 


0 


0 


/o 


0 


m 1 




-e 2 


-e 3 


0 


0 


0 


0 


fo. 


_m 8 



where e 0 — b^f^ — b 2l ,e l —b l2 ,fi — b 22 ,e 2 ~^ > wfi ~ ^23^3 ~ ^i4»/*3 — ^24' 
and for (31) we obtain 



e 0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


-e 4 


0 


e 3 


-c 2 


/l 


m l 


0 




0 


0 


0 


0 


0 


0 


0 


0 


e 4 


0 


-«3 


0 




-/ 2 


m 2 


0 


0 


e 0 


0 


0 


0 


0 


0 


0 


e A 


0 


0 






0 


/a 


m 3 


0 


0 


0 




0 


0 


0 


0 


e 4 


0 


0 


0 


/1 


"/ 2 


"/a 


0 


m 4 


0 


0 


0 


0 


e 0 


0 


0 


0 


0 


e 3 


-e 2 


-*I 


0 


0 


0 


-/ 4 


m 5 


0 


0 


0 


0 


0 


e 0 


0 


0 




0 


fx 


-/ 2 


0 


0 


/ 4 


0 


m 6 


0 


0 


0 


0 


0 


0 


e 0 


0 




-fx 


0 


"/a 


0 


-f* 


0 


0 


m-, 


0 


0 


0 


0 


0 


0 


0 


e 0 




fi 


/a 


0 


/ 4 


0 


0 


0 


m 8 


0 


0 


0 


-f* 


0 


"/a 


fl 


fx 


/o 


0 


0 


0 


0 


0 


0 


0 


m 9 


0 


0 


-h 


0 


"/a 


0 




~e 2 


0 


fo 


0 


0 


0 


0 


0 


0 


m 10 


0 


-/ 4 


0 


0 


/ 2 




0 


-e 3 


0 


0 


/o 


0 


0 


0 


0 


0 


m n 


h 


0 


0 


0 


/1 


^2 




0 


0 


0 


0 


/o 


0 


0 


0 


0 


m 12 


0 


/a 


fi 


-ex 


0 


0 


0 


-e A 


0 


0 


0 


0 


/o 


0 


0 


0 


m 13 


"/a 


0 


f 


^2 


0 


0 


e 4 


0 


0 


0 


0 


0 


0 


fo 


0 


0 


m 14 




/l 


0 


e 3 


0 


-e 4 


0 


0 


0 


0 


0 


0 


0 


0 


fo 


0 


m 15 




^2 


-e 3 


0 


e 4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


fo _ 


m X6 



= 0 



(38) 

where e 0 — b n ,f 0 — b 2v e l — b l2 ,f x —b 22 ,e 2 —b l3 ,f 2 —b 23 ,e 3 —b u ,f 3 —b 24 ,e 4 —b l5 ,f 4 —b 25 ; 
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Similarly, systems of linear equations with larger sizes could be obtained for the equation (24-3) 
(where s = 1,2,3,... , r = 2). 

In general, the size of the quadratic matrices of these matrix equations (corresponding with the second 
order equations of the type (24-3), i.e. for r = 2 ) is X x 2 s . But exceptionally, this size is reducible to 
2 s " 1 x 2 s ~ l (only for the case of the second order equations), as it was for equations (29) - (31). 

Generally, the size of the quadratic matrices of the matrix equations (corresponding with the general cases 
of the equation (24-3)) is r s x r s . For all values of s , r , these matrix equations (according to (23), and 
corresponding to the general cases of the equation (24-3)) are derivable and calculable. 



Meanwhile, as we previously noted, any of quadratic matrices (i.e. the matrices At in (23)) in equations 
(32) - (38) and so on, is just one of the possible matrices, from which we may obtain. Definitely, there are 
other matrices (different ones, but with the same size and similar structure) which we may obtain for 
constructing other equivalent cases (based on (23)) of the equations (31) - (38); we just selected those 
individual matrices (in (32) - (38)) because of their particular structures, as well as their particular 
applications that will be in the next section. 



As some examples of the third order cases of equation (24-3), the systems of linear equations 
corresponding with two equations 

zri* M =*iAAi. (39) 

2 3 

= b nKK + b l2 b 22 b 22 ; (40) 

q=\ p=\ 

respectively, are 





0 


0" 


m 1 






0 


fo 


0 


m 2 


= 0, 


(41) 


0 


0 


So_ 


m 3 







10 



0 


0 


0 


e l 


0 


0 


0 


0 


e o 


m l 


0 


0 


0 


0 


/, 


0 


fo 


0 


0 


m 2 


0 


0 


0 


0 


0 


Si 


0 


So 


0 


m 3 


0 


0 


e o 


0 


0 


0 


e l 


0 


0 


m 4 


fo 


0 


0 


0 


0 


0 


0 


/l 


0 


m 5 


0 


So 


0 


0 


0 


0 


0 


0 


Si 


m 6 


«1 


0 


0 


0 


0 




0 


0 


0 


m 1 


0 


/. 


0 


fo 


0 


0 


0 


0 


0 


m 8 


0 


0 


Si 


0 


So 


0 


0 


0 


0 


m 9 



where e Q =b n ,f Q =b 2V g Q =b 3V e l =b n ,f =b 22 ,g l =b 32 . 



The size of the quadratic matrix in the matrix equation, corresponding to the next 3 r order equation, i.e. 

q=\ p=l 



Concerning the fourth order cases of equation (24-3), such as 



Sn^ 9 = ^iAAAi> (43) 

q=\ p=l 

2 4 

ZIl^ = b n b iAAi + b 12 b 22 b 32 b 42 ; (44) 

q=\ p=l 



respectively, the systems of linear equations corresponding to them are 



e 0 


0 


0 


0" 


m { 


0 


fo 


0 


0 


m 2 


0 


0 


So 


0 


m 3 


0 


0 


0 


\_ 


_m 4 
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0 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


e 0 


m l 


0 


0 


0 


0 


0 


fl 


0 


0 


0 


0 


0 


0 


fo 


0 


0 


0 


m 2 


0 


0 


0 


0 


0 


0 


Si 


0 


0 


0 


0 


0 


0 


So 


0 


0 


m 3 


0 


0 


0 


0 


0 


0 


0 


h 


0 


0 


0 


0 


0 


0 


K 


0 


m 4 


0 


0 


0 




0 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


m 5 


fo 


0 


0 


0 


0 


0 


0 


0 


0 


fl 


0 


0 


0 


0 


0 


0 


m 6 


0 


So 


0 


0 


0 


0 


0 


0 


0 


0 


Si 


0 


0 


0 


0 


0 


m 7 


0 


0 


K 


0 


0 


0 


0 


0 


0 


0 


0 


K 


0 


0 


0 


0 


m 8 


0 


0 


0 


0 


0 


0 


0 


e 0 


0 


0 


0 


0 




0 


0 


0 


m 9 


0 


0 


0 


0 


fo 


0 


0 


0 


0 


0 


0 


0 


0 


/l 


0 


0 


m 10 


0 


0 


0 


0 


0 


So 


0 


0 


0 


0 


0 


0 


0 


0 


Si 


0 


m n 


0 


0 


0 


0 


0 


0 


K 


0 


0 


0 


0 


0 


0 


0 


0 


K 


m 12 


-«l 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


e 0 


0 


0 


0 


0 


m 13 


0 


/l 


0 


0 


0 


0 


0 


0 


fo 


0 


0 


0 


0 


0 


0 


0 


m 14 


0 


0 


Si 


0 


0 


0 


0 


0 


0 


So 


0 


0 


0 


0 


0 


0 


m 15 


0 


0 


0 


*. 


0 


0 


0 


0 


0 


0 


K 


0 


0 


0 


0 


0 


_ m i6 


e e Q 


= b ivfo 




So 


b 3V h 0 




b AV e x 


= b n>fi = 


b 


Si = 


Z? 32 , /ij 


= b A2 . 









= 0 



(46) 



Similarly, as we noted above, for the fifth and the higher order cases of (24-3), with the larger number of 
unknown elements, we get the matrix equations including the quadratic matrices with the size r s xr s . 

Meanwhile, we can use the following linear relations (as the general rules) for transforming the second, 
the third, the forth and the higher order cases of the form (24) into the form (25); e.g. for the second order 

z%rk=znv W7> 

ii ,1*2=1 p=l q=\ p=l 

we have 

s s s 

b ll = C V K =2X2^2' ^12 = C 2> ^22 = E 5 2/ 2 C , 2 ' — K= C s' b 2s=T, B sh C i2> < 48 ) 

/ 2 =1 ; 2 =1 ; 2 =1 

and for the third order: 
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i % 3 rR=£ri*v (49) 

h,l 2 ,l 3 =l p=l q=l p=l 



we have 



b u c p b 21 Cj, b 3l ^ I B Ui c i , b l2 c l , b 22 c 2 , b 23 » 

/ 3 =i / 3 =i 



^=c P ^=c„ & 3 ,= jx^,.- V-. + i) =c ^ V^d =c i' *Wi)=Z B *i* c . 

/ 3 =1 i 3 =1 

s 



! ' 3=1 (50) 



Similarly, for transforming the forth order and the higher order cases of equation (24) into (25), we can 
define some linear transformations such as (48) and (50). However, if necessary, it is possible for 
transforming (24) into (25), we define other linear transformations, as well. 

According to the particular applications of the second order cases of the equation (24-3) in the next 
section (the section 3.), here we do consider and introduce some of the properties of the matrix equations 

(34), (36), (37), (38). 

First, we consider the following equation (which its left part is a special case of (26)): 



Y j B ij {c i c j -d i d j ) = Q (51) 



i,j=0 



Where B tj - B fi . Now using the matrix equations (34) (for ft - 1), (36) (for ft - 2), (37) (for ft - 3 ) and 

n 

(38) (for ft = 4) and so on, as well as the linear transformations of the type e { = /B ii jc i +dj), 

;'=o 

ft =c i -d i ,ot 
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f 

Jo 






A 




c, —d x 


f 3 




^2 ~ d 2 


Jn_ 




c — d 

_ n n_ 



(52-1) 







B QQ 


B 0l 


B 02 ■ 


■ B 0n 


c 0 + d Q 


e l 




B 1Q 


B n 


B 12 ■ 


■ K 


C, + J[ 






B 2Q 


B 2l 


B 22 ■ 


■ B 2n 


c 2 + d 2 






_ B nO 




B n2 ■ 


B nn _ 


C n +d n_ 



(52-2) 



obtain the following matrix equations that correspond to equation (51), respectively 



[5 00 (c 0 +J 0 )][m 1 ] = 0, 

i 

Y, B oj( c j +d j) c i~ d 



7=0 



m, 



= 0, 



^Bojicj+dj) 

7=0 

0 



0 



7=0 
2 



7=0 7=0 



7=0 7=0 



c 0 d 0 
0 



(53) 









mj 


(c 2 -d 2 ) 


m 2 


0 


m 3 




m 4 


c -J 
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0 


0 


0 


0 


-e 3 


e 2 


/r 


mj 


0 


^0 


0 


0 


^3 


0 


-«i 


/ 2 


m 2 


0 


0 




0 


"«2 


«i 


0 


/ 3 


m 3 


0 


0 


0 


e o 


-/l 




"/a 


0 


m 4 


0 


-/a 


/ 2 


*1 


/o 


0 


0 


0 


m 5 


f 3 


0 


-/l 


^2 


0 


/o 


0 


0 


m 6 


-fl 


/l 


0 




0 


0 


/o 


0 


m 7 


-«1 


"«2 


-e 3 


0 


0 


0 


0 


/o. 


_m 8 



= 0, 



(56) 



where 



e o =ll B oj( c j+ d j)> fo= c o- d o> 

7=0 

3 

7=0 
3 

e 2=^ B 2j( C j +d j)> f 2 =C 2 -d 2 , 
7=0 

3 

e 3 =Z 5 3 i ( c J f 3 =C 3 ~d 3 . 



(56-1) 



7=0 
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0 


0 


0 


0 


0 


e o 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


e o 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


-/ 4 


0 


0 


0 




0 


-f 


0 


-h 


0 


0 


u 


h 


0 


0 


0 


fl 


0 


fs 


fl 


-e. 


0 


-u 


0 


f 


^2 


0 


fl 


fx 


0 


^3 


0 


-ex 


?1 


-e 3 


0 


e A 



where 



0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


e 3 


" e 2 


~ e l 




0 


0 


e 3 


0 


/, 


-fl 


0 




0 


~ e i 


-A 


0 


-f 


0 


0 


e o 




f 2 


fs 


0 


-f 3 


fl 


/, 


fo 


0 


0 


0 


0 


e l 




0 


fo 


0 


0 


~ e x 


0 




0 


0 


fo 


0 




e 3 


0 


0 


0 


0 


fo 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


-e 4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




IX 


(Cj+dj), f Q =c Q 


-d 0 , 




e x = 


4 

2X 

7=0 


(Cj+dj), 


fl = c l - 


-d v 




e 2 = 


7=0 


(Cj+dj), 


fl = c l 


— d 2 , 




e 3 = 


7=0 


(Cj+dj), f 3 =c 3 


— d 3 , 




e A = 




(Cj+dj), f A =c A 


-d A . 





0 


e 3 


~ e l 


/, 


m l 


-e 3 


0 


~ e l 


-fl 


m i 


~ e i 


~ e i 


0 


fs 


m 3 


fx 


~fi 


-f 


0 


m 4 


0 


0 


0 


-/ 4 


m 5 


0 


0 


/ 4 


0 


m 6 


0 


-/ 4 


0 


0 


m 1 


n 


0 


0 


0 


m 8 


0 


0 


0 


0 


m 9 


0 


0 


0 


0 


m 10 


0 


0 


0 


0 


m n 


0 


0 


0 


0 


m xi 


fo 


0 


0 


0 


m 13 


0 


/o 


0 


0 


m 14 


0 


0 


fo 


0 




0 


0 


0 


fo _ 


_ m !6 



= 0 



(57) 



(57-1) 



We must note that there are not the same linear transformations such as (51-1) - (51-2) (as exceptionally, 
they exist for (51)), for the third and the higher order equations of the form 

n n 

Z B vk ( c i c i c k - d l d j d k ) = 0 , Z B uki ( c i c j c k c i ~ d,d ] d k d l ) = 0 , . . . . (58) 

ij,k=0 i,j,k,l=0 

Now, by the following choices, 
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B = 



B oo 


B Q1 


B Q2 . 


■ ■ B 0n 


B 10 


B n 


B 12 ■ 


■ ■ K 


B 2Q 


B 2\ 


B 22 ■ 


■ ■ B 2n 



nO nl '-'nl 



C = 



c o 




d Q 








fo 






d l 








A 


C 2 




d 2 








fs 




, D = 




, E = 




, F = 








d n 








fn 



(59) 



we can rewrite the transformations (52-1) and (52-2) as follows 



E = B(C + D), F = C-D, 



(60) 



from (60) we also get 



C = Y 2 (B- l E + F), D = y 2 (B l E-F). 



(61) 



Where B is the inverse of B ; according to equation (51), B is a symmetric matrix; also we assume 
that det5 7^0. The relations (60) and (61) indicate that there is an one to one correspondence between 
the components of C, D and the components of E, F ; using this property, as well as the relations (60) 
and (61), we will determine the solutions of the system of linear equations (54) - (57), on the basis of the 
solutions of equations (34), (36), (37) and (38). 

Utilizing the standard and specific methods of solving the systems of linear equations in the set of 
integers [7], respectively, we obtain the following sets of the solutions for equations (34), (36), (37) and 

(38), (for unknowns e { and / ( ). 

Firstly, we must emphasize that the general and standard solution of the equation of the type 
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1=1 



n-1 



=0 (62) 

in the set of integers, is as follows [7, 8]: 

x t =a n k i , (i = 1,2,3,..., n-1), x n = ^a i k i (63) 

i=i 

where the parameters k i are arbitrary integers, and where we supposed " a n ^ 0 ". 

Now, for the equations (34) we get the following solutions (where we supposed " m 2 ^ 0 "): 

e Q -k 2 m 2 , f^-k^, e l -k 1 m 2 , f x --k 2 m x (64) 

where the parameters k v k 2 ,m l ,m 2 are arbitrary integers. 

For the equations (36) we have (where we supposed " m 4 ^0 "): 



(65) 



where the parameters k l ,k 2 ,k 3 ;m l ,m 2 ,m 3 , m 4 are arbitrary integers. 

Using Remark 2., we can also get the following general solutions for the system of linear equations (36) 
(where we supposed " m 4 ^ 0, k 4 ^ 0 "): 

e Q = k 3 m 4 - k 4 m 3 , f Q = k 2 m x -k x m 2 , e, = k 2 m 4 - k 4 m 2 , 

(66) 

/j =k l m 3 - k 3 m x , e 2 = A;,m 4 - fc 4 mj , / 2 = k 3 m 2 - k 2 m 3 . 
where the parameters k l ,k 2 ,k 3 ,k 4 ;m l ,m 2 ,m 3 ,m 4 are arbitrary integers. 

For the equations (37), the following solutions are obtained (where we supposed " m g ^ 0 ", and there is 
also a condition for the parameters m i (see below)): 
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(67) 

where the parameters k x ,k 2 ,k 3 ,k 4 are arbitrary integers and the parameters m ; should satisfy the 
following equation (as a necessary condition for the parameters m ; that exist in the parametric solutions 
(67)): 

m 4 m 8 = -mjm 5 - m 2 m 6 - m 3 m 7 (68) 

Since the parameter m 4 does not exist in the solutions (67), the condition (68), easily, could be solved by 
the following choices: 

w 8 =1, m 4 = -U\U 5 - u 2 u 6 - u 3 u 7 , 

m l =U\, m 2 =u 2 , m 3 =u 3 , (69) 

m 5 =u 5 , m 6 =u 6 , m 7 =u 7 . 

Now using Remark 2. and the relations (69) and (67), the solutions of (37) are determined as follows 
6q — k^t , ^0 — k^ti^ ~\~ ^2^2 ^1^3 ? — k-^t } — k^n^ ~\~ k^ti^ k^U-j * 

(70) 

^2 — kryt ^ — k^i&2 ^1^5 k^u^i , — k]t j — k^ti^ kr^xi^ k^ti^- 

where the parameters t,k l ,k 2 ,k 3 ,k 4 ;u l ,u 2 ,u 3 ,u 5 ,u 6 ,u 1 are arbitrary integers, and (t 0) . 

Similarly, we obtain the following solutions for the matrix equation (38) (where we supposed " m 16 ^ 0 ", 
and there are also five conditions for the parameters m i (see below)): 



e Q - k 5 m l(> , f 0 - k 4 m y - k 3 m 2 + k 2 m 3 - k x m 5 , e { - k 4 m l6 , f y - -& 5 m, + k y m l2 - k 2 m l4 + k 3 m l5 , 
e 2 =k 3 m l6 , f 2 =k 5 m 2 +k l m n -k 2 m i3 -k 4 m l5 , e 3 =k 2 m l6 , f 3 = -k 5 m 3 -k { m {f) + k 3 m u + k 4 m l4 , (71) 
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where the parameters k v k 2 ,k 3 ,k 4 ,k 5 are arbitrary integers and the parameters m ; should satisfy the 

following equations (as the necessary conditions for the parameters m i that exist in the parametric 
solutions (71)): 

m 4 m 16 = m[m 13 + m 2 m 14 - m 3 m 15 , 
m 6 m 16 = m l m ll + m 2 m i2 -m 5 m 15 , 

m 7 m i6 =_m i m io _m 3 m i2 + m 5 m i4' ( 72 ) 
m 8 m 16 = -m 2 m lQ + m 3 m n -m 5 m l3 , 

m 9 m i6 =-m io m i5 + m ii m i4 — m 12 m 13 . 

In like manner, since the parameters m 4 ,m 6 ,m 7 ,m & ,m g do not exist in the solutions (71), the conditions 
(72), will be solved by the following choices: 

m l6 = 1, 

TYL^ — W^Wj j "I - W 2^12 ^5^15' 
TYtq — ^1^10 W-^Wj2 W^Wj 
TYLq — Z^Wjq "I - M-^Wjj ^5^13' 

(73) 

YY\/q — Z-/ j 0^ 1 5 ^ 1 1^ 14 ^ 1 2^ 13' 

m x =u x , m 2 =u 2 , m 3 =u 3 , 

^5 — ^5' "^10 — ^10' "^11 — ^11' 

m 15 = w 15 . 

Using the relations (71) and (73), as well as Remark 2., the solutions of (38) become as follows 

£q — ^5^' — k^A^ k^U^ k]H^ i — k^tj — k^iy k]H^2 ^2^14 ^3^15' 

— k-^tj — ^5^2 k^U'Y y ^2^1 3 ^4^15 5 ^3 — k^t •> — ^"5^3 ^1^10 ^3^1 3 k^tt^ ^ * (^~J 1 - 1 ^ 

— k-^t j — ^5^5 ^2^10 ^3^11 ^4^12* 

where the parameters t,k v k 2 ,k 3 ,k 4 ,k 5 ;u v u 2 ,u 3 ,u 5 ,u lQ ,u lv u l2 ,u l3 ,u l4 ,u l5 are arbitrary integers, and 
(f*0). 
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Similarly, the parametric solution of the matrix equation (with the size 32x32) corresponding to 
equation, 

!>,/,. =0, (74) 

;=o 

similar to (64), (65), (70), (71-1), will be gotten, but with sixteen additional conditions for parameters m i , 

(these conditions include sixteen homogenous second order equations, that each equation contains only 
four terms, similar to (68) and (72)). These conditions could be solved, easily, with some specific choices 
for parameters m i , similar to (69) and (73). In general, the parametric solution of the matrix equation (i.e. 
the system of linear equations) corresponding to the second order equation of the form 

n(n + 1) 

will lead to (2" 1) number of conditions for parameters m i (including the four terms 

homogenous second order equations), that these conditions will be solved by some specific choices for 
parameters m i , similar to the above choices (the choices (69) and (73)), and ultimately, the general 
solution of (75) will be obtained. 

Meanwhile, the solutions (64), (65), (70), (71-1) could also be represented as follows, respectively 







k 2 t 






e 1 




k x t 






7o" 




' 0 




V 


X 




- W] 


0 





(76) 



where we just supposed ( m 2 = t) and ( m x =u l ). 















e l 




k 2 t 












_k l2 t_ 








/o" 




" 0 




-u 2 




/l 






0 


u 3 


k 2 


fl_ 




u 2 


-u 3 


0 





(77) 



where we just supposed ( m 4 = t ) and ( m x = u x , m 2 -u 2 , m 3 =u 3 ). 
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e 0 




k 4 t 
















k 3 t 
















k 2 t 


•> 














k x t 












/o" 




o 






u 2 


u 3 


"*4 


fi 




- u x 


0 


- w 7 


u 6 


^3 


h 




- u 2 


u 1 


0 


-u s 




u. 




_-u 3 


-u 6 


W 5 


0 





(78) 







fe 5 ? 












































e 3 
































/o" 




" 0 




- u 2 


u 3 


""5 




/l 




-Mj 0 


u l5 


-u u 


M 12 


^4 


/ 2 




u 2 




0 


- u l3 


M n 


^3 


/ 3 




U 3 M 14 




0 


- w 10 


^2 


A. 




_ u 5 


M 12 


-u u 


M 10 


0 


_*. 



(79) 



Corollary 1. If we define the matrix K as 



K = 



(80) 



where we suppose K 0 , then using (60) and (61) we get the following sets of the relations 
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C = / 2 (tB- l +U)K, 
D = / 2 (tB 1 -U)K, 

(81) 

E = tK, F = UK, 
K*0, detfi^O; 

D = (tB l -U)(tB l +U) l C, 

E = 2t(tB 1 +uy 1 c, 

F = 2M(tB 1 + Uy l C, , (82) 
K = 2(tB 1 +uy 1 c, 
C*0, detfl^O; 

C = (tB 1 +U)(tB 1 -U) l D, 
E = 2t(tB 1 -U) l D, 

F = 2M(tB 1 -U) l D, (83) 
K = 2(tB l -U) l D, 
D*0, detfl^O. 

In point of fact, the formulas (81) are the (integer) parametric solution of equation (51), where the 
matrices B, C, D, K,U have been defined by the relations (59) and (76) - (80), and the parameter t is an 
arbitrary integer (t ^ 0). On the other hand, the formulas (82) and (83) show that, if we suppose that the 
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values c i (or d i ), (where i = 1,2,3,..., n) are given values, then we can calculate the values d, (or c,. ) 
in terms of them and the matrices B,U . 

Remark 3. The condition (68) (of the parameters m i ) concerning the parametric solutions (67), and the 

conditions (72) (of the parameters m i ) concerning the parametric solutions (71), could be solved by other 
methods as well. 

Since the parameter m 4 does not exist in the solutions (67), the condition (68), easily will be solved by 
the following choices (for the parameters m i ): 

m 4 =0, (84-1) 

m 8 : a free Integer parameter (m 8 ^ 0), (84-2) 
m l m 5 + m 2 m 6 + m 3 m 7 = 0. (84-3) 

where equation (84-3), according to the solutions (65) and (66), also could be solved as follows 
(including two sorts of the solutions): 

m l = u 3 v 4 , m 2 = u 2 v 4 , m 3 = u 1 v 4 , 

m 5 = u 2 v l -u 1 v 2 , m 6 =u 1 v 3 -u 3 v v m-, = u 3 v 2 -u 2 v 3 ; (85) 
m 4 = 0, m 8 : a free Integer parameter (m 8 ^ 0); 

and 

m l = u 3 v A - u 4 v 3 , m 2 =u 2 v 4 -u 4 v 2 , m 3 =u l v 4 -u 4 v l , 

m 5 =u 2 v l -u 1 v 2 , m 6 =u l v 3 -u 3 v l , m 7 = u 3 v 2 -u 2 v 3 , m 4 =0, (86) 

m 8 : a free Integer parameter (m 8 ^ 0). 
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where the parameters U V U 2 ,U 3 ,U 4 ,V V V 2 ,V 3 ,V 4 are arbitrary integers. By replacing the values of m j , 

(from the relations (85) or (86)) in formulas (67), as well as taking into account formulas (84-1) and (84- 
2), we get a new general parametric solution for the matrix equation (37). 

As another similar case, since the parameters m 4 ,m 6 ,m 1 ,m & ,m 9 do not exist in the parametric solutions 

(71), the conditions (72) (for the parameters m i , existed in the solutions (71) of the matrix equation (38)) 
could be solved by the following choices, as well 



m 4 =m 6 =m 7 =m s =m 9 = 0, (87-1) 

m 16 : a free Integer parameter (m 16 ^ 0), (87-2) 

m 1 m lQ + m 3 m l2 -m 5 m 14 = 0, (88-1) 

m 1 m l j + m 2 m l2 - m 5 m l5 = 0, (88-2) 

m 1 m l3 + m 2 m l4 -m 3 m 15 =0, (88-3) 

m 2 m lQ + m 5 m l3 -m 3 m n =Q, (88-4) 

m 1( /n l5 + m l2 m l3 -m ll m l4 = 0. (88-5) 



As the equation (88-5) could be derived from (88-1), (88-2), (88-3) and (88-4), we will not consider it in 
the next relevant calculations. 

Referring to the solutions (65) and (66) (for equation (36) that corresponds to the quadratic equation 
(29)), respectively, the equations (88-1), (88-2), (88-3) and (88-4) are solved as follows, 

first, using the solutions (65) we get: 

m 1 = u 4 v 5 , m 13 = w 3 v 2 -u 2 v 3 , 

(88-1) —> m 2 =u 3 v 5 , m u = u 2 v 4 -u 4 v 2 , (89-1) 

-m 3 =u 2 v 5 , m 15 =u 4 v 3 -u 3 v 4 ; 
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m l = u 4 v 5 , m 1 j = u 3 v l - WjV 3 , 
(88-2) —> m 2 =u 3 v 5 , m l2 =u l v 4 -u 4 v l , (89-2) 

-m 5 = Ul v 5 , m l5 =u 4 v 3 -u 3 v 4 ; 



m x =u 4 v^ m l0 = (-u 2 )v l -u l v' 2 , 
(88-3) —> m 3 =-u 2 v 5 , m l2 =u l v 4 -u 4 v l , (89-3) 

-m 5 =u { v 5 , m l4 =u 4 v' 2 -(-u 2 )v 4 ; 



m 2 =u 3 v 5 , m 1Q = u 2 v[-(-u l )v 2 , 
(88-4) —> -m 3 =u 2 v 5 , m n -(-u l )v 3 -u 3 v[, (89-4) 

m 5 =-u,v 5 , m l3 = u 3 v 2 -u 2 v 3 . 

By the definitions of the type v[ = —v v v' 2 = —v 2 , the solutions (89-1) - (89-4) could be simplified, for 
instance. 

Now using the solutions (66), we get the other type (more expanded one than (89-1) - (89-4)) of the 
solutions for equations (88-1) - (884), respectively, 

m l = u 4 v 5 -u 5 v 4 , m l3 -u 3 v 2 -u 2 v 3 , 
(88-1) ->• m 2 =u 3 v 5 -u 5 v 3 , m 14 = u 2 v 4 -u 4 v 2 , (90) 

-m 3 -u 2 v 5 -u 5 v 2 , m l5 =u 4 v 3 -u 3 v 4 ; 
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m l - u 4 v 5 -u 5 v 4 , m ll = u 3 v l —u l v 3 , 



(88-2) -> 



m 2 =u 3 v 5 -M5V3, m l2 = u l v 4 -u 4 v l , 



(91) 



■m 5 =M 1 v 5 -M 5 v 1 , m 15 =u 4 v 3 -u 3 v 4 ; 



m, = w 4 v 5 -w 5 v 4 , mj 0 = (-m 2 )v 1 -m 1 (-v 2 ), 



(88-3) -> 



m 3 = -(u 2 v 5 -w 5 v 2 ), m l2 = u 1 v 4 -u 4 v 1 , 



(92) 



-m 5 =1/^5 -m 5 v p m 14 = w 4 (-v 2 )-(-w 2 )v 4 ; 



m 2 = w 3 v 5 -w 5 v 3 , m 10 =m 2 (-v 1 )-(-Mj)v 2 , 
(88-4) -» - m 3 = u 2 v 5 - u 5 v 2 , m l , = (-w, ) v 3 - m 3 (-Vj ), (93) 

m 5 = -(u 1 v 5 -MjVj), m 13 = m 3 v 2 -u 2 v v 

By simplifying the relations (90) - (93), ultimately we get the following set of the general solutions for 
equations (88-1) -(88-4): 
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m 1 = u A v 5 - u 5 v A , m 2 = u 3 v 5 -u 5 v 3 , 

m 3 =u 5 v 2 -u 2 v 5 , m 4 =0, 

m 5 = u 5 v [ - u l v 5 , m 6 = 0, 

m 7 =0, m 8 = 0, m 9 = 0, 

m l0 = u l v 2 -u 2 v 1 , m n =u 3 v l -u l v 3 , 

m l2 = u l v 4 -u 4 v l , m l3 = u 3 v 2 -u 2 v 3 , 

m lA = u 2 v A -u A v 2 , m l5 =u A v 3 -u 3 v A , 

m 16 : a free Integer parameter (m 16 ^ 0). (94.) 



where the parameters U 1 ,U 2 ,U 3 ,U 4 ,U 5 ;v l ,V 2 ,V 3 ,V 4 ,V 5 are arbitrary integers. By replacing the values of 

m i (from the relations (94), in formulas (71)), as well as taking into account formulas (87-1) and (87-2), 
we get another form of the general parametric solution for the matrix equation (38). 

Similarly, for the systems of linear equations corresponding to (75), with more variable elements (i.e. 
larger values of n in (75)), the similar conditions and relations to (84-1) - (84-3) and (87-1) - (87-2) and 
(88-1) - (88-5) and so on, could be chosen, ultimately. Then, based on them, we can determine the 
solution for the system of linear equations corresponding to each specific case of the quadratic equation 
(75). 



3. Deriving the General (and Unique) Structure of the Laws Governing the 
Fundamental Forces of Nature. 



If we assume that in the special relativity condition, the components of the n -momentum are discrete, i.e. 
have integer values in the invariant and the energy-momentum relations 



g"P,P v =g M P,P' v , (95) 
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uv / \2 00/ m Q C \2 / 0/r s 

8 P M P v =(-m 0 c) = g i-j=) (96) 

where g^are some constant coefficients, and p„,p' are the components of the momentum vector in 

two reference frames, then the relations (95) and (96) are the special cases of the relation (51). Hence, 
using the matrix relations (53) - (57), we get the following systems, respectively, where they correspond 

to (95) and (96) ( S i are the integer parameters similar to the parameters m i ), 
for (95) we have: 

k°>o + K)]h] = 0 (97) 



8 (Pv+Pv) Pl~ Pi 

-g 1V (Pv + Pl) Po-Po_ 



= 0 



(98) 



where V — 0,1 ; 



8 0V (P V + Pl) o 

o g 0v ( Pv + Pl) 
Pi -p'l Pi- p'l 

-g lv {p v +p' v ) g 2v (p v +p'J 



g 2v (p v +p' v ) 
g iv (p v +p' v ) 

Po-Po 
0 



Pi- Pi 
(Pl-p'l) 

0 

Po-Po 



= 0 



where V = 0,1,2 ; 
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(100) 



(100-1) 
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e 0 = 8 0l, (Pv + P'yl fo = Po-Po> 

ei=g lv (P v + P' v ), fi = Pi-p[, 

e 2 = g 2v (P v + Pi), f 2 = P 2 -p' 2 , 

e 3 =g 3v (p v + P' V X f 3 = P 3 -Pr 



(100-2) 
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where we have 



S 4 S \6 — ^1*1 3 ~*~ "^2^14 ^lS' 



^6^16 — S X S X1 + ^2^12 ^5^15' 



S 7 S 16 — ^lO ^3^12 ^5 ^14' 



(101-1) 
(101-2) 
(101-3) 
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S & S 16~ S 2 S 10 + S 3 S l 1 S 5 S 13> 



S 9 S 16 ~ S \Q S \5 +S \l S U S l2 S iy> 



(101-4) 
(101-5) 



e 0 =8 V (Pv + P'vX fo = Po-p'o, 



ex=g' v (P v + P' v ), fi = Pi-p[, 



e 2 = S v (Pv + pI), f 2 = P 2 -p' 2 , 



(101-6) 



^ =g V (P V + P' V \ fi=P3-P^ 



e 4 = g V (P V + P' V ), f 4 = P 4 -p' 4 . 



m n c 



And for (96) we get (where we suppose jU^O: p' M = 0, p' 0 = /=^)' respectively 
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(102) 



00/ m o c 



Pi 



= 0 



(103) 



where V = 0,1 ; 
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0 
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-g lv p v 
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p v -g°\- 

Px 



g 2v p v 



~g 2v p v 



g Xv P v 
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. m Q c . 
+ ("/==) 



= 0 (104) 



where V = 0,1,2 ; 
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= 0 (105) 



where we have 

s 4 s s + s l s 5 + s 2 s 6 + s 3 s 1 = 0 (105-1) 



e =s Qv n -g°°( m ° C ) f =p +( m ° C ) 



e x =g v p v , f = Px, 

(105-2) 
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(106) 

where we have 



S A S 16 = S 1 S 13 + S 2 S u -S 3 S l5 , (106-1) 

Vl6=Vll+Vl2-Vl5. ( 106 " 2 ) 

s 7 s 16 = -j^,, - s 3 s 12 + 5 5 5 14 , (106-3) 

■^le = — ^2^10 + — s 5 s w (106-4) 

^9*1 6 = ^10^15 ^1 1^14 ^12^13' (106-5) 
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Qv 

^0=8 Pv 

e 2 =g 2v P v > f 2 = P2> (106-6) 
e 4 =g 4v p v , f 4 = P 4 - 

Remark 4. In this section we will use the geometrized unites [9], Einstein notation, as well as the 
following conventions: 

-Metric sign convention: (+,-,-,...,-), 

r p =8 r p + r p r z -d r p -r p r z 

-Riemann and Ricci tensors: 

R = —R v 

-Einstein tensor: (R /lv -y 2 Rg /lv ) = -MT /lv + ... . (107) 

Corollary 2. It is worth to note that concerning the relations (97) - (101) for the components p and p' M 

, we can use the relations (81) - (83) and the solutions (76) - (79), to determine the linear transformations 
between two reference frames. In other word, the general forms of the linear transformations between two 
reference frames, directly and immediately, are determined from the relations (97) - (101) for all 
dimensions. 

Now, we use equations (102) - (106) for deriving the equations governing the fundamental forces of 
physics. 

First, we consider the following definitions for the canonical operators: 

p M ^>ihV M (108) 
g f,v (const.) — > g MV (components of the metric tensor) (109) 

s ; . — > s t =F MV , Z , ... (components of the field tensor) (110) 



.00/ m 0 C \ f _ _ , 00/ m 0 C 



g G=)> fo = Po+g m '(- 



fi=Pv 
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We do substitute the operators (108) - (109) with the corresponding quantities in equations (102) - (106). 
According to the structures, the compositions and the number of equations of the systems (102) - (106), 
we do conclude that there exist only three kinds of anti-symmetric tensors that their components could be 

substituted with the parameters S t (except in the equation (102), which is a special and trivial case), and 

transform the systems (102) - (106) to the tensor equations. These tensors are a 2 nd order tensor, a 3 rd 
order tensor and a 4 th order tensor. The 4 th order tensor of these tensors just matches the Riemann tensor 

Rpvpa (which at the same time, it is necessary for calculating and specifying the components of the 
metric tensor g MV ), the other two tensors could be written as Z and F , that should be anti- 
symmetric with respect to the indices u, V such as: Z = — Z , F ' — —F . 



Thus firstly, corresponding to equations (102) - (104) (that they don't contain any condition for 
parameters s i ), we get the following tensor equations, respectively 

(the tensor equation corresponding to relation (102) is a special and trivial case, that Riemann tensor also 
vanished there; in any case we just write it down, assuming a tensor such as F substituting with 5 X and 

where we assume g 00 = 1) 

(102) -> D* M F M =0 (Hl-1) 



where jU = 0, and g 00 = 1, s 1 —>s 1 =F 0 . 



D [p F MV] =0, (112-1) 

(103)- 

K FM v=-Jl E) ( 112 - 2 ) 



where p,ju,v- 0,1, and s l ->s 1 =F 10 , s 2 ^s 2 = <p (E> , J ( V E) = -D v <p (E) . 



*VW=0, (112-3) 

(103)- 

D;Z" vp =-J^ (112-4) 

where p, a,p,V = 0,1 , and ^ -> ^ = Z lQp , s 2 -> s 2 = (p p N) , J™ = -D v (p p N) . 
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*VW=0, (H2-5) 

(103)- 



where A, A cr, //, V = 0,1 , and s x -> ^ = /? 10yC(T , s 2 -> s 2 = p£> , = -£> v ^ ) 



D [p F MV] =0, (113-1) 

(104)- 

DlF" v =-J™ (113-2) 

where p,jU,V - 0,1,2 , and 

-> = F 10 , s 2 -> s 2 = F 02 , s 3 ^s 3 = F 2V s 4 ^s 4 = cp (E \ J ( V E) = -D v (p {E) . 



*VW=0, (113-3) 

(104)- 

D*fi%=-J%> ( 113 - 4 ) 

where p, a,/J,V — 0,1,2 , and 

s l —> 5] = Z 10p , s 2 — » s 2 = Z 02 ^, s 3 — > s 3 = Z 2lp , s 4 — > s 4 = ) = —D v <p p 



*VW= 0 ' ( 113 - 5 ) 

(104)- 

^V=-C (H3-6) 

where /I, /?, <J,/U,V - 0,1,2 , and 
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s^s^R^, s 2 ->s 2 =R 02pa , s 3 ^s 3 =R 21fXT , s 4 -*s 4 =(p£\ J™=-D v q>™ 



where in equations (111-1) - (113-6) we have 



D M =V M + ^K, (114-1) 



h M 



D- M =V M -?±k M . (114-2) 



ju = 0: k 



n 

l 



ju*0: k = 0, 



V v /(^ =0 , /fWf'-^y^, (116-1) 



(115-1) 



v/r=Q' c=c-^^ z v (H6-2) 



v v / ( ?=o, C=C-yW- ^v ( ?=o. (116-3) 



In principle, we suppose and conclude that the unique tensor equations (112-1) - (113-6) (including the 
definitions and relations (114-1) - (116-3)) are the laws governing the fundamental forces of nature, 
respectively in the dimensions D — 2,3 . 



In the system of linear equations (102) - (104), the parameters S t were just arbitrary integer parameters, 
that we substituted the components of tensors F^ v ,Z Mvp ,R MVfXT with them (based on the principal 
operator definitions (108) - (110)). 

Before writing the tensor equations corresponding to the relations (105) and (106), which are similar to 
the tensor equations (112-1) - (113-6) (as we will show), let at first mention the following remark and 
corollary. 



Remark 5. For the next systems of linear equations, i.e. (105) and (106) (and so on), the situation for the 
parameters S { is a bit different. There are some conditions for the parameters S. (including the quadratic 
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equations (105-1) and (106-1) - (106-5)), that should be considered and solved. In the previous section, 
we dealt with this situation in two different ways. We had two types of the general solutions for these 
conditions; one includes the solutions of the form (69) and (73) (that simply, they respectively become the 

solutions of the conditions (105-1) and (106-1) - (106-5), by m i — > s i ). The other one includes the 

solutions of the forms (85) and (86) (that simply, they become the solutions of the condition (105-1), by 

m i — >s i ), and the solutions of the forms (89-1) - (89-4) and (94) (that simply, they become the 

solutions of the conditions (106-1) - (106-5), by m i — > s t ). In all these solutions for the conditions 

(105-1) and (106-1) - (106-5), the parameters s i , necessarily, are written in terms of the new parameters 

U i and v { , which we may represent them as the general form 



where the parameters U i and V ( . are arbitrary integers. Now, for dealing with this situation, in principle, 
we accept that: 

" Due to the tensor components (corresponding to the operators s t , according to the definitions (108) - 
(110)) substitute with the parameters S { , and at the same time the parameters S t are not arbitrary 
parameters and they defined by (117) (where u j and V ( . are arbitrary parameters), then we conclude that 
the operators s t should have the following form, as well 



where U i and v i are some operators that substitute with the parameters U t and V, , and they should be 

specified for the tensor components corresponding to the operators S t . So, from the above conditions and 
arguments, in principle, we accept and conclude that 

[s, = H i (u v u 2 ,u 3 ,...,u m ;v v v 2 ,v 3 ,...,v n )]^[s i =^,.(M 1 ,M 2 ,M3,...,M m ;v 1 ,v 2 ,v 3 ,...,v n )]." (119) 

Meanwhile, according to the parametric solutions (69), (73), (85), (86), (89-1) - (89-4) and (94), the form 
H t is not a unique form. In Corollary 3. by using (119), the form H i , as well as U i and V, could be 

specified (in fact, beforehand, we will clarify that which form(s) of H '. are acceptable), such that they 



s t =H i (u v u 2 ,u 3 ,...,u m ;v v v 2 ,v 3 ,...,v n ) 



(117) 



s t =H i (u v u 2 ,u 3 ,...,u m ;v 1 ,v 2 ,v 3 ,...,v n ) 



(118) 
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will be consistent with the operators S t (that correspond to the components of the tensors 

p 1 7 J? \ 

fiV //vp' nvpo>' 

Corollary 3. Concerning Remark 5. and the relations (118) and (119), for specifying the form H t , and 
the operators U i and V, in the structures of the tensors F /Jv ,Z /Jvp ,R flvpa ., basically, we will use the 
Riemann tensor as the base, as on one hand it is necessary for specifying the components of the metric 
tensor g M v , and on other hand it is a basal tensor with a "specific structure". 

Thus, on this basis, concerning the form H ) in (117), (118) and (119), only the parametric solutions of 
the type (86) (with respect to m i —> s t ): 

s 1 =w 3 v 4 -M4V3, s 2 = u 2 v 4 -u 4 v 2 , 

S 3 = WjV 4 - M 4 Vj , S 5 = U 2 V l - WjV 2 , 

(117-1) 

S 6 = U l V 3 ~ U 3 V l ' S 7 =U 3 V 2 -U 2 V 3 , S 4 =0, 

5 8 : a free Integer parameter (j g ^ 0). 
for equation (105-1), as well as only the parametric solutions of the type (94) (with respect to m i — > s t ): 

S l =M 4 V 5 -M 5 V 4 , S 2 =U 3 V 5 ~U 5 V 3 , 

s 3 = u 5 v 2 -u 2 v 5 , s 4 =0, 
s 5 =u 5 v l -u l v 5 , s 6 =0, 
s 7 =0, s s = 0, s 9 = 0, 
s l0 =u l v 2 -u 2 v l , s n = u 3 v l -u l v 3 , 
s l2 =u l v 4 -u 4 v l , s l3 = u 3 v 2 -u 2 v 3 , 
s u =u 2 v 4 -u 4 v 2 , s 15 =u 4 v 3 -u 3 v 4 , 

s 16 : a free Integer parameter (s l6 * 0). (117-2) 
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for equation (106-1) - (106-5), are acceptable; furthermore, by starting from the Riemann tensor and its 
specific structure in (107) and using the relation Y x = g px T p , we get 



Rpvv = (^r w , -ry^ ) - (d M v pav ) (120) 

If we define the operator C M , operating on a second order tensor Y vp as follows 

C M K P = a(d v Y pM + d p Y VM - d M Y vp ) (121) 
where a is a non-zero constant, then we may rewrite (120) such as 

R wv = y a K d v C P -^*A)«*. -(V, -r^cj^] (122) 

Now, we define the operator D as well 

D pM =(d M C p -T A pM C z ) (123) 
Then the relation (122) for the Riemann tensor could be written as follows 



R = 1/(D s -D s ) (1241 

papv /Q pvocjfi ppoov' V ' 



Meanwhile, the natural generalization of the operator C , operating on an arbitrary n th order tensor 



According to (122) and (124), the only and the most general form representing the structure of the 
Riemann tensor, is 

R M vpa = R pap v = C^/nAp ~ ^pfPov ) ( 1 26 ) 
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where A =D , B =8 . (Ill) 



Now, as stated in Corollary 3., and Remark 5., we use and apply the relation (126) (as a basal and 
principal criterion), as well as the relations (117) - (119), for determining the general formulation of the 

operators 5, (that substitute with the parameters s i in the systems (103) - (106)), which correspond to 

the components of the Riemann tensor. That is, respectively 



(103): {(as l =aR 10piT = A p0 B ITl -A pl B IT0 , $ 2 =<pfJ\ 

(as 1 = A 0 B 1 -A 1 B 0 , (128) 
s 2 : a free Integer parameter (s 2 ^ 0)) } ; 



(104) : { (as, = aR lQpa = A pQ B al - A pl B a0 , as 2 = aR 02pa = A p2 B a0 - A p0 B a2 , 



as 3 = aR 2Xpa = A pX B a2 - A p2 B aV s 4 = <p p G J ), (as, = A Q B, - A X B Q , 



s 4 : a free Integer parameter (s 4 * 0)) } ; 



(129) 



(105) : { (as, = aR lQpa = A p0 B al - A pl B aQ , as 2 = aR 2Qpa = A p0 B a2 - A p2 B a0 , 

as 3 = aR 3Qpc7 = A pQ B a3 - A p3 B a0 , s 4 = 0, as 5 = aR 23pc7 = A p3 B a2 — A p2 B a3 , 

as 6 = aR 3lpy = A pl B a3 -A p3 B aV as 7 = aR l2fXT = A p2 B al -A pl B a2 , s s = p£>), 

(as, = AqB, - A,B Q , as 2 = A Q B 2 - A 2 B Q , as 3 = A Q B 3 - A 3 B 0 , s 4 = 0, 

as 5 = A 3 B 2 - A 2 B 3 , as 6 = A X B 3 - A 3 B, , as-, = A 2 B X - A,B 2 , 
s & : a free Integer parameter (s s ^ 0)) } ; 
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(130) 



(106) : { (as 1 = aR l0pc7 = A pQ B al - A pl B aQ , as 2 = aR Q2/X7 = A p2 B a0 - A pQ B a2 , 

a h = aR 3Qpa = A pQ B a3 - A p3 B aQ , s 4 = 0, as 5 = aR Q4pcj = A p4 B aQ - A p0 B a4 , 
s 6 = 0, 5 7 = 0, s & = 0 s g = 0, as 1Q = aR 34/X7 = A p4 B a3 - A p3 B a4 , 
^ ! = aR 42pa = A p2 B a4 - A p4 B a2 , as 12 = aR 4Xpa = A pl B a4 - A p4 B al , 
as 13 = aR 23pcj = A p3 B a2 - A p2 B a3 , as l4 = aR 13/XT = A p3 B al - A pl B a3 , 
as l5 = aR 2lpa = A pl B a2 - A p2 B al , s 16 = (p% y ), 

s 4 = 0, as 5 = A 4 B Q - A 0 B 4 , s 6 = 0, s 7 = 0, s s = 0, s g = 0, 
as 1Q = A 4 B 3 - A 3 B 4 , as n = A 2 B 4 - A 4 B 2 , as l2 = A X B 4 - A 4 B { , 

(131) 

s 16 : a /ree Integer parameter (s l6 ^0))}. 

where the parameters A. and B are arbitrary. In the above formulas for S { , the element a is just an 
arbitrary parameter ( a ^ 0 ). But in the above formulas representing S t , the element a will be specified 
later (that is a — iti ). 

As we stated in Corollary 3., as well as according to formulas (128) - (131), it is clear that for the form 
H i in (117), (118) and (119), only the relations (117-1) and (117-2) are acceptable, which are consistent 
with the formulas (126) and (128) - (131). Firstly, it is easy to show that the algebraic formulas 
representing s { in (128) and (129), are also consistent with previous algebraic conditions of the 

parameters 5 ; in the systems (103) and (104) (where they just were arbitrary integers). For example, 

regarding the algebraic formulas representing S { in (129), by the following choices 
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(132) 

where the parameters , s[ are arbitrary integers, directly we can show that the parameters S,. are also 
arbitrary and each could take any integer value. 

Secondly, concerning the formulas (130) and (131), directly we could show that by the following choices, 
respectively, 

A 0 =av 4 , B 0 =u 4 , A l =av 3 , B x =u 3 , 

A 2 =av 2 , B 2 =u 2 , A 3 =av 1 , B 3 =u 1 , (133) 
s 4 - 0, 5 8 : a free Integer parameter (s s * 0); 



Aq = av 5 , B Q = u 5 , Aj = av 4 , B x =u A , A 2 = -av 3 , 
B 2 =-u 3 , A 3 =-av 2 , B 3 = -u 2 , A 4 =av x , B A =u v 

(134) 

5 4 =0, 5 6 = 0, 5 7 =0, 5 8 = 0, 5 9 = 0, 

5 16 : a /ree Integer parameter (s l6 *0) 



the formulas (130) and (131) completely accord with formulas (117-1) and (117-2) (and only with them, 
as it should be). 

Now, according to the above suppositions and approaches, such as Remark 5. and Corollary 3. and the 
relations (130) and (133), we derive the following unique field equations as the laws governing the 
fundamental forces of nature in the fourth dimension, respectively 

D [p F MV] =0, (135-1) 

(105)- 
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where p,jU,V- 0,1,2,3 and 

5j — > S 1 = F lQ , S 2 — > 5 2 = ^20' ^3 — ^ ^3 = ^30» 

s 4 — > s 4 = 0, s 5 — > s 5 = F 23 , s 6 — > s 6 — F 3l , 



D lp Z flvia =0, (135-3) 

(105)- 

Dp%,=-J™ (135-4) 



where p,(T,jU,V — 0,1,2,3 and 

5j — > iSj = Z l0p , s 2 — > 5 2 = Z 20p , 5 3 — > s 3 = Z 30 ^, 

^4 ~ ^ ^4 = 0, ^5 ~ > ^5 = Z 23p , 5 6 — > S 6 = Z 3 lp , 



Sl ^s,=Z np , s^s,=(p p N \ J^=-D v cp^. 



*VW=0, (135-5) 

(105)- 



Where A, p, <7,jU,V = 0,1,2,3 and 
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s 4 ^>s 4 =0, s 5 ^s 5 =R 23fXT , s 6 ^s 6 =R 31/XT , 
s _>s = # c _>s = G / G > J (G) =-D 0 {G) 



Where the definitions and relations (114-1) - (116-3) are applied to the equations (135-1) - (135-6) as 
well 1 . 

Corollary 4. Now, as stated in Corollary 3. and Remark 5., according to the relations (117) - (119) and 
(123), (125), (127) and (128) - (131), as well as (126) (as a basal and principal criterions for the 

specifying the structures of the quantities S. , that correspond to the components of the field tensors), the 
following relations regarding the structures of the field tensors F (the representer of the 
electromagnetic field) and Z (the representer of the strong nuclear field, as it should be) in all 
equations (112-1) - (113-6) and (135-1) - (135-6), are necessary, respectively 

F MV = -(D MV Q-D Vfl Q) (136) 
a 

Z Mvp = Sw V ^iP^H^-D^H,) (137) 

where Q and H are a scalar field and a vector field, and where we have a = i%, as well as the 
following gauge operators that ultimately, are obtained for the fields F v and (or S 0l v ) : 



F V +— — A 



A M = C M Q = -d M Q, 



(138-1) 



is (E) 

F MV = V V A M ~ VA +JL ^l\'\l d38-2) 



( 

where g is the coupling constant; 



1. Similarly, using the relations (131), (134) and (108) - (110), we may obtain the tensor equations corresponding to the system 
(106). However, there might be appeared some limitations for the higher dimensions if we suppose the discreteness of the other quantities such 
as space-time coordinates and so on. 
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and 

V^V^^— H M , 

(139 " 1} 

L MV = V v H M -V M H v+ ^-[H v ,H M ], 
where g <A ° is the coupling constant. 



Corollary 5. According to the unique structures of equations (135-1) - (135-4), we can conclude that 
magnetic monopoles do not exist in nature. 

Remark 6. The Generalized form of Einstein field equations could be obtained from equations (1 12-6), 
(113-6), (135-6). Using the (second) Bianchi identity and the conventions (107) we get 



V R" =V R -V R (140) 

fi vpcr p va a vp V > 



Now from (140) and equations (112-6), (113-6), (135-6) and the assumption 

J vpa =-%K{VJ vp -V p T va ) + %7tB{V a Tg vp -V p Tg va ) (141) 
where 7^ v is the stress-energy tensor (T = T M M ), and g MV is the metric tensor and B is a constant, we get 



lift 

R MV = -^v" BT S MV ) -~^ K »v- 4S, V (142) 



where q is a constant (that appears naturally, when we obtain equation (142)), and K — V M k v (where 

k v has been defined in (115-1)), and K — —K Vfj , K M M = 0 (due to its complex coefficient); Equation 
(142) for (112-6) (two dimensional case), specially, takes the following form 
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^T MV + -^K MV = / 2 (SxT)g MV , (143-1) 
R MV = - l A(SxT)g MV + y 2 A gfJv ; (143-2) 



R-A = -SttT (144) 



where A-2q, 5 = 0 and A is the cosmological constant. 

Regarding equation (113-6) (three dimensional case), the specific form of the field equation (142) is as 
follow 

itn 

R MV ~ YiRg, v = -%xT MV - ~yK MV - Ag fiv (145) 
where A = / 2 q, 5 = 1 and A is the cosmological constant. 

And concerning equation (135-6) (four dimensional case), the field equation (142) takes the following 
specific form, as well 

ifti 

R MV -y2Rg MV = SxT MV —^K MV -Ag MV (146) 

where A = q, B = j/^ and A is the cosmological constant. 

In equations (142) - (146) we may have the (total) conservation law as follows 1 

V(T MV + -^Z flv ) = 0 (147) 

" h 



1. Meanwhile, according to the second Bianchi identity, from the equations (113-5) and (135-5) we may have 
Ijfl 

£1, b ^ 0 '. R-abrxr = ^ ' " lat m ^ e case ^0 ^ ^ ^ as mvariant mass °f a supposed gravitational force carrier particle), there may 

fl 

appear some additional conditions for the components of the Riemann tensor. 
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Remark 7. In the special relativity conditions, equations (112-1) - (112-4), (113-1) - (113-4) and (135-1) 
- (135-4) also could be written as follows 



(ifia fJ d fJ -Im 0 )[F] = 0, 
(itia fJ d fJ -Im Q )[Z] = 0. 



(148) 
(149) 



where / is the identity matrix; and where for equations (112-1) - (1 12-4) (i.e. two dimensional case) we 
have 



[F] = 



"1 0 " 




' 0 f 




, a = 




0 -1 




-1 0 



F 



10 

(E) 



, [Z] = 



(150) 



and for equations (113-1) - (1 13-4) (i.e. three dimensional case) we have 



[F] 



1 0 

0 1 

0 0 

0 0 

F 

1 no 



Fr 



21 

(£) 



0 
0 

-1 

0 



0 
0 
0 



a 



[Z] = 



7 



0 0 0 1 
0 0 10 
0-100 
-10 0 0 



a 



0 0 10 
0 0 0 -1 
-10 0 0 
0 10 0 



(151) 



and for equations (135-1) - (135-4) (i.e. four dimensional case) we have 
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0 


0 
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0 


0 
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1 
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0 
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0 
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0 


-1 


0 


0 


0 


0 


0 


1 


0 




0 


0 


0 


„ 1 
, a = 


0 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


1 
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0 




-1 


0 


0 
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0 
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0 


0 
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0 




0 


-1 


0 




0 


-1 


0 


0 
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0 


0 


0 


0 


0 


0 


0 


0 




0 


0 


-1 




-1 


0 


0 


0 


0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


-1 


cf 




" 0 


0 


0 


0 


0 


1 


0 


0" 


0 


0 
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0 


0 
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0 


1 
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0 


-1 
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0 
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0 
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0 
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-1 
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0 
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-1 
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0 


0 
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1 
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0 
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0 




-1 
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1 


0 




0 
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0 


0 
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0 
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0 


0 


1 


0 


0 


0 


0 


0 


-1 




0 


0 


0 


0 


0 


0 




0 


0 


-1 


0 


0 


0 


0 


0 



[F] 







7 


F 

1 20 




7 


F 

1 30 




7 


0 

F 

1 23 


, [F] = 


0 

7 

^23p 








Fn 




7 

^I2p 


q> (E) 







(152) 



All these matrices (including previous matrices) correspond to Clifford algebras [10]. 

Corollary 5. In some special cases such as: (m 0 ^ 0) and (<p (E> - 0, <p =0), the equations (148) and 

(149) (as well as the equations (112-1) - (113-6) and (135-1) - (135-6)) are turned into some (wave) 
equations similar to Pauli and Dirac equations. Which we can suppose that they are the corrected and 
generalized forms of Pauli equation (in two dimensional case), and Dirac equation (in three dimensional 
case) and so on. 
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4. Conclusion 



In the section 2., since the set of algebraic axioms (17) - (23) for integers, have been formulated in terms 
of the quadratic (n x n) matrices (with arbitrary "n"), we may conclude that for complete representing 
the algebraic property of integers, necessarily and sufficiently, the quadratic matrices (nxn): 
(K]„x„»[^y]„x„»[ c «>]nx„v e Z nxn ) should be applied. The (old) algebraic axioms (10) - (16-2) that had 
been formulated in terms of the elements: a x ,a v a^... € Z, (that in fact, we can assume that these 

elements are lxl matrices such as: [a] M (= a), [b] Ui (= b), [c] M (= c),... e Z M (= Z)) are not sufficient 
for complete introducing the algebraic properties of integers. 

In the section 3. of this article we derive the general and unique forms of all field equations governing the 
fundamental forces of nature, including equations (112-1) - (113-6) and (135-1) - (135-6). This 
derivation was based on a new mathematical approach (represented in the section 2., concerning the 
algebraic properties of the domain of integers), as well as the assumption of discreteness of the 
components of the relativistic n-momentum. The obtained laws include three separate types of the field 
tensor equations (and only three), representing the force of gravity, the electromagnetic force, and the 
(strong) nuclear force (for dimensions D >2). In some special conditions, these equations could be turned 
into the (wave) equations that are similar to Pauli and Dirac equations. In actual fact, they are the 
corrected and generalized forms of the current field equations including Maxwell equations, Yang-Mils 
equations and Einstein equations, as well as (in some special conditions) Pauli equation, Dirac equation, 
and so on. In addition, according to obtained laws, we also showed that magnetic monopoles could not 
exist in nature. The results obtained in the section 3. demonstrate the efficiency of the theory of 
linearization (and its algebra) and a wide range of its possible applications. The proposed mathematical 
structure in the section 2., doubtless can be useful in many fields where the discreteness of the relevant 
quantities are supposed and applied. 
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